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Holographic and other imaging systems utilizing coherent light introduce 
a speckled or noise-like pattern in the image of a diffuse object which 
severely degrades image quality. It is desirable to understand this effect 
quantitatively. Intelligent design in many cases requires knowledge of 
the mean-square value, spatial power spectral density, and autocorrelation 
junction of the noise-like fluctuations. These quantities have been deter- 
mined for the image of a uniform diffuse object. Major results are: 

(i) The mean-square value of the fluctuation in the image intensity 
is equal to the square of the mean intensity. 

(ii) One can decrease the relative magnitude of the noise-like fluctua- 
tions at the cost of a corresponding increase in the aperture required of 
the optical system (or hologram) over that required to resolve the desired 
image in a spatial frequency sense. In a holographic facsimile or TV 
system, this calls for a corresponding increase in electrical bandwidth. 

(Hi) The improvement in (ii) is not possible for direct viewing with 
the human eye, since the resolution of a healthy eye is known to be limited 
by diffraction at the iris. 

I. INTRODUCTION' 

Holographic and other imaging systems using coherent light have 
been receiving considerable attention lately. 1 ' -• 8 - 4 Most analyses on 
this subject assume that the object reflects specularly, or transmits 
specularly if the object is a transparency, i.e., the reflectivity or 
transmissivity of the object varies smoothly. Most objects, however, 
are more nearly diffuse reflectors. When the image of a diffusely re- 
flecting object is formed it will be covered with a noise or grain-like 
structure 5, "• 7 which is the speckle pattern which one sees when laser 
light is used to illuminate an object. 
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In this paper we investigate the noise-like or speckled nature of 
the image of a uniform diffuse surface. It should be emphasized that 
we are interested in the properties of the image in contradistinction 
to the direct backscattered field studied by Goldfisher. 8 We show 
that the intensity consists of two parts. The first is the mean or en- 
semble average intensity and is proportional to the intensity which 
would be obtained if incoherent light were used for illumination. This 
is the desired component of the image and might be likened to a 
signal. The second part of the image is the speckled or noise-like 
component which tends to obscure the average intensity. This noise- 
like component occurs because of the random phase angles associated 
with the scattering centers comprising the microstructure of the dif- 
fuse surface. The spatial autocorrelation function and power spectral 
density of the speckle pattern in the image arc found, and are shown 
to be dependent upon the size of the aperture stop. It is shown that 
the variance of the intensity fluctuation is equal to the square of the 
mean intensity. The fluctuation may be reduced, however, if one is 
willing to sacrifice resolution by recording the image on film whose 
resolution is much poorer than that set by the aperture of the optics. 
Unfortunately, this alternative is not available when viewing with the 
human eye, since the resolution of a healthy eye is known to be de- 
termined by the diffraction limit of the iris." This seems to place 
definite limitations upon the use of coherent light in visual systems. 

II. ARBITRARY APERTURE 

The model which we shall use for a diffuse object is shown in 
Fig. 1. Although the object is shown to be a granular transparency, 
it could equally well have been shown as a reflector without loss of 
generality. The essential point is that a monochromatic coherent light 
wave of unit intensity is assumed to be scattered by a random set of 
point scattercrs. Each scatterer is assumed to be a unit scatterer which 
is many wavelengths in depth from its neighbor. The relative phase 
of the wave scattered from each scatterer may be assumed to be a 
random variable which is statistically independent of the phase of the 
waves scattered from other scattercrs. Any phase change between 
and 2tt is equally likely. Multiple scattering will be neglected. 

The scattered field just to the right of the granular transparency 
can be expressed by the equation 

Fo(x, y) = E Kx -x it y- ,/,y 9 '", (1) 
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Fig. 1 — A uniform wave of coherent light is incident on a transparency com- 
posed of randomly distributed unit point scattercrs. Light collected by the 
aperture A, placed in the far-field, is imaged by lens L on plane P. 

where #,- is the relative phase of the wave scattered from the scatterer 
located at x = x i} y = ?/,-. 8 h .r, and iji are assumed to be random vari- 
ables uniformly distributed in the intervals (0,2tt), (— X,+X) and 
( — Y,+ Y), respectively. Notice that because of our assumptions, the 
statistics of the scattered field are independent of any deterministic 
variation in the phase of the illuminating field. 

A Fourier transform relationship exists between the scattered field 
given by (1) and its far-field. The far-field is given by 



Fi(Z, v) = f f " F (x, yy 

J — OD J _-o 



(2i/X./)(j-£ + i/t|) 



dx dy 



(2) 



__ y it i + j(2»/Xd)(xi{ + »,i)) 



where we have suppressed the time factor c +,! ". Notice that each scat- 
terer has produced a plane wave, and that the slope of the phase front 
of each wave with respect to the $, 77 axes is determined by the position 
(Xi , yi) of the random scatterer. 

Let the far-field /*\(£, -q) be passed through an aperture having an 
amplitude transmission function //(£, 17), and then through a lens which 
is placed a distance 2 behind the aperture. Since the field at the back 
focal plane of a lens is a Fourier transform like function of the field in 
front of the lens, an image of the granular transparency, as modified 
by the aperture, will be formed in the back focal plane, and is given by 10 
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where c = tt(z - f)/\f, and where h{t, u) and #(£, v) are Fourier 
transform pairs in the sense 

Kt, U) - f° f " fffo f y»-««*^ # r/rj. (3b) 

J —oo J— ao 

Notice that except for the unimportant phase factor e ,i:< ' } , (3a) 
differs from (1) for the field at the granular transparency only in that 
h( ) functions have replaced the delta functions. That is to say, the 
delta function of light from the scatterer at (x { , y t ) is reproduced as a 
broadened h( ) function located at v = — (f/d)x, ; , co = — (//d)y.- • 
The image is reversed, and magnified by the factor m = f/d. Notice 
that because of the random phase 0,- of each, the impulse functions will 
add vectorially in a random fashion when they overlap one another. 

The situation is analogous to passing shot noise impulses through 
a low-pass filter having an impulse response h( ). The impulses are 
broadened into h{ ) pulses whose width depends inversely upon the 
filter bandwidth. In the coherent light case, however, the process is 
two dimensional and the applied impulses have random phase angles 
distributed uniformly between and 2-rr, rather than being constrained 
to be positive impulse functions as is the case for shot noise. 

The quantity of greatest interest to us is the intensity of the image, 
which is found by multiplying the image field by its conjugate. 

I{v, co) = F 2 (v, u>)F%( V , co) (4) 

~ h h w + xrf ' x/ + u) ll \\j + xd ' x/ + \dr 

The uniform diffuse object is assumed to exist in the region — X ^ 
x ^ +X, — Y ^= y ^ +Y. The number K of point scattcrers in this 
region is a random variable, as are their positions (x,, y t ) and their 
relative phase angles ( . We may, therefore, obtain the ensemble 
average of the image intensity / by averaging (4) with respect to the 
2/v + 1 random variables consisting of the K positions {x,, ;/;>, K 
phase angles 0,-, and K itself: 

£* *L i Mi. 

vX/ ' Xrf ' X/ ~ l ~ \d 



'-£"-£"[5*fe- 
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?/._-> 



V h*(— 4- -^i — -4- 

frt \X/ "■" Xd ' X/ T Xrfy 

W(&i i V\ ; ^ , y 3 ; • • • a* > 2/* ; 0i ; • • • ^ ; K) 

dxx dijt •■■ dx K dy K dd x ; • • • dd K dK, (5) 
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where W( ) is the multi-dimensional probability density function. 
Now the positions (x,, y,) are considered to be statistically inde- 
pendent variables, as are the relative phase angles 6,. They are also 
independent of K, so we may simplify (5) to obtain 

7 - 
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(6) 

We see that the above expression vanishes unless 0, = 6 k , i.e., i = k. 
Further, all of the h( ) functions have the same shape so that if the 
size of a resolution element in the image is small compared to the 
field of view, i.e., the extent of h(v/\f, io/\f) is small compared to X 
and Y, then we may replace the limits of integration ±X and ±F by 
± oo to obtain 

/.^jTmioa ( 7) 

p\ (u, v) is the autocorrelation function of the aperture impulse 
function h($, 77) , i.e., 



Pi(w, v) = j J h*((, r)h{t + u, t + w) r// r/r 



(S) 



If we now assume that the number of scatterers per unit area of the 
transparency has a Poisson distribution of mean A", then the mean 
intensity is 

7 = ^A\Yp,(0, 0). (9) 

Next we wish to determine an expression for the autocorrelation 
function of the intensity, from which we may determine the spatial 
power spectral density and variance of the noise-like fluctuations. The 
autocorrelation function of the intensity as given by (4) is 



R(r, = l 2 (p, a>)I 2 (p + r, « + t) 



(10) 
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Because of the statistical independence of the phase angles ft, posi- 
tions (xi, yi) and K, and because of the assumed uniform distribution, 
we may simplify (10) to 

R(t, t) = f"w(K)dK 

' hi h hi h i-oc 2X J-oo 2F ■ " J_ M 2X Loo 2 7 
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We see that the integral vanishes unless 

i = k and n = m 
or 

n = k and i = m, 
which gives 



(11) 
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Now, we have two subcases here. There are K(K—1) terms for which 
k ¥= m, and there are K terms for which k = m. 

R(t, = / " K(K - l)W(K) dK 
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(13) 



Assuming that the distribution of scattcrers W(K) is Poisson and 
using the definition of h ( ) given in (3b), straightforward evalua- 
tion of the integrals in (13) yields 



R(t, t) = T 



1 + 



pJj/jK t/jX) I 2 ' 
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where />i (a, v) is defined in (8) and 
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p,(u, V) = r° f " | h(r, I) | 2 | h{r + tt, / + u) 



p 2 (r//X, f//X), (14) 



</t </« 



= autocorrelation function of the magnitude squared of 
the aperture impulse function. 

The spatial power spectral density is found by taking the Fourier 
transform of (14). After simplification we obtain 



S(q,p) = j** jT" R(t, &-»•<«**» d, 



dt 



(i\f 



H(Xjq, X/p) 



I 2 \8(<i,p) + ^ L j\H(Xf< h X(p) | 2 
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where © stands for convolution. In particular we define 
F(\fq, \/p) ® G(\fq, X/p) 

- ( + ° [ +X F*(x, y)G(\fq + x, X/p + y) dy dx. (16) 

J — oo "— oo 

Equation (9), which gives the mean intensity of the image, and 
(15), which gives the power spectral density of the intensity fluctua- 
tions, are the major results of this section. 

III. CIRCULAR APERTURE 

Now consider the special case of a circular aperture of radius r n 
and let it be located on axis so that 

where 

r = + V? + i? 2 . 
The average intensity in the image plane is given by (9) and is 

7 = d 2 \*N Pl (0, 0) = T.\(\dr r f, (18) 

where pi (0,0) was evaluated from the integral 

Pi(0, 0) = f^ r | Hfo v) I 2 di d V = ttt 2 . • (19) 

J — oo •'—00 

Evaluation of the integrals in (15) gives the power spectral density 



S(q, V) - I' 



1 L 2 . _, / s \ 2 / s \ ( s 



*•»+*{*-;** is) -;\sNi-\s: 



(20) 



1 /, 2 . _, / s \ 2/s| / _^ 

2s f , 

where 

g, p = image plane spatial frequencies in rectangular coordinates, 

s = + v7 + p 2 

Sr = r c //X = cutoff frequency produced by diffraction at the circular 
aperture. 

2wr- c 
= overlap factor. 
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The overlap factor F warrants some discussion. Basically it is equal 
to the average number of point scatterer image centers contained 
within an area equal to that occupied by the image of a single 
scatterer. That is, a single point scatterer located at (0,0) in the 
object plane would produce a point in the image plane at (0,0) 
having intensity 



= h 



!-/i[fey7+j' 



IT V " + M " 



The intensity is down 11 approximately 50 percent at (2irr e /f\) 
y/v\ + u> 2 = "\/2, and the area covered by the image of the point 
scatterer at this 50 percent value is A = ir{v\ + u 2 ) = f\ 2 /2irr 2 . For 
a diffuse object, the average number of imaged scattering centers per 
unit area in the image plane is n = (d/j)'X. If we define the overlap 
factor F as the average number of scatterer image centers falling in the 
area of one of these images we have 



F = nA = 



d 2 \\\ 



For a truly diffuse surface, the overlap factor F » 1 so that (20) re- 
duces In 



8(q, v) = T' 



.^rt+St 1 --™^ 



■K \ZS, 
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(21) 



which is plotted in Fig. 2. Note that it is symmetrical about the 
vertical axis. For very small spatial frequencies, (21) can be approxi- 
mated bv 



S(q.p) = / 2 [5(f / ,p)+^]. 



(22) 



The total fluctuation or noise power occurring in spatial frequencies 
less than some frequency Si is 



'' ^ /T '- 



(23) 
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IV. CONCLUSIONS 

We have found that the image of a uniform diffuse object illumi- 
nated with monochromatic coherent light consists of two parts. The 
first is the mean or ensemble average given by (18), for a circular 
aperture, and is proportional to the intensity which would be obtained 
if noncoherent light were used for illumination. This is the desired 
component and might be likened to the signal component of the 
image. The square of this term appears as the first term in (20), 
(21), and (22), and as the delta function in Fig. 2. The second part 
of the image is a grainy or noise-like component which tends to 
obscure the mean intensity or signal. This noise-like component 
occurs because of the random phase angles associated with the point 
scatterers comprising the microstructure of the diffuse object. This 
component is shown as the second term in (20), (21), and (22), and 
as the continuous part of the power spectrum in Fig. 2. Integration of 
(21) shows that the variance of the noise-like fluctuations in the 
intensity is equal to the square of the mean intensity (or to the signal 
power). This is fortunate to the extent that when the signal is small, 
the noise is likewise small. However, while our analysis was for the 
particular case of a uniform diffuse surface, we can safely predict 
that for nonuniform diffuse objects fine detail in the image will be 
largely obscured by the noise-like fluctuations if resolution is limited 
by diffraction. 

The noise-like fluctuations in the image can be reduced if one 
records the image on film whose modulation transfer function has a 
bandwidth which is much smaller than the diffraction limit of the 
optical system. The high-frequency noise in Fig. 2 will not be resolved 
in this case. For instance, if one requires the "signal-to-noise" ratio 
to be increased from unity to 10 3 (30 dB), then from (23) we see that 




!M=«i*^-*--'(il-J-i(^)fU? 



(S/2S 



Fig. 2 — Section of the spatial power spectral density for a uniform diffuse 
surface imaged through a circular aperture. The complete two dimensional 
spectrum is obtained by rotating the above curve about the vertical axis. 
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the diffraction bandwidth s,. of the improved optical system must be 
10 3/2 = 31.6 times the bandwidth S| resolvable by the film, and there- 
fore by the whole system. (Since most transducers produce a signal 
which is proportional to the intensity of the incident light, it seems 
appropriate to consider the square of the mean intensity as signal 
power and the variance of the intensity fluctuations as noise power.) 

Although we have analyzed the very special optical system shown 
in Fig. 1, our results are not critically dependent upon the placement 
of the aperture. The aperture could be the lens aperture rather than 
an independent physical device, or it could be the aperture defined 
by the finite size of a hologram, for instance. Our results should also 
hold approximately for the human eye, since the resolution of a 
healthy eye is known to be determined by the diffraction limit of the 
iris. 9 The predicted value of unity for the signal-to-noise ratio is the 
right order of magnitude for what one observes when laser light is 
used for illumination if one is careful to hold the eye stationary and 
hence not average the noise out as a function of time. Although 
moving the eye tends to average out the noise, the residual noisiness 
remains objectionable. This places definite limitations upon the use 
of coherent light in visual systems. 

The author wishes to thank Messrs. C. B. Rubinstein and A. B. 
Larsen for helpful discussions. 
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